Law of corresponding states for osmotic swelling of vesicles 
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As solute molecules permeate into a vesicle due to a concentration difference across its membrane, 
the vesicle swells through osmosis. The swelling can be divided into two stages: (a) an "ironing" 
stage, where the volume-to-area ratio of the vesicle increases without a significant change in its 
area; (b) a stretching stage, where the vesicle grows while remaining essentially spherical, until it 
ruptures. We show that the crossover between these two stages can be represented as a broadened 
continuous phase transition. Consequently, the swelling curves for different vesicles and different 
permeating solutes can be rescaled into a single, theoretically predicted, universal curve. Such 
a data collapse is demonstrated for giant unilamellar POPC vesicles, osmotically swollen due to 
the permeation of urea, glycerol, or ethylene glycol. We thereby gain a sensitive measurement of 
the solutes' membrane permeability coefficients, finding a concentration-independent coefficient for 
urea, while those of glycerol and ethylene glycol are found to increase with solute concentration. In 
addition, we use the width of the transition, as extracted from the data collapse, to infer the number 
of independent bending modes that affect the thermodynamics of the vesicle in the transition region. 



I. INTRODUCTION 

Membrane vesicles are made of a closed bilayer of am- 
phiphilic molecules in aqueous solution, having length 
scales of 0.1-100 /im. Used as simplified models of bio- 
logical membranes, they have been one of the most ex- 
tensively studied systems in soft- matter physics [l|, 
Vesicles usually enclose both solvent (water) and solute 
molecules. Such vesicular capsules are ubiquitous in cell 
functions [3| and used as microreactors and delivery vehi- 
cles in various biomedical and cosmetic applications 

The hydrophobic core of the bilayer membrane poses a 
kinetic barrier to the permeation of water and water- 
soluble molecules into and out of the vesicle. Conse- 
quently, the permeability coefficients of various solutes 
across various membranes span a very wide range of val- 
ues @ . For example, the permeability coefficients of wa- 
ter, urea, and potassium ions through the membranes 
used in the current work are of the order, respectively, of 
10 2 , 10~ 2 , and 10~ 8 /xm/s [5|]. This implies that the time 
scales for permeation of water and the solutes examined 
here are well separated. Thus, the vesicles can safely be 
assumed to reside in a semi-permeable regime, in which 
their volume quickly (essentially immediately) re-adjusts 
through osmosis to a change in the number of enclosed 
solute molecules. As more solute molecules enter, there- 
fore, the vesicle progressively swells — first approaching 
a spherical shape (the "ironing" stage), subsequently in- 
flating as a sphere (the stretching stage) , and eventually 
rupturing (osmotic lysis) @. 

In an earlier theoretical work we argued that a vesicle 
in such a semi-permeable regime should reach the end of 
the ironing stage (i.e., the maximum volume-to-area ra- 
tio) critically, through a continuous transition j?:]. That 
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theory was restricted to unstretchable membranes. Here 
we extend the theory to the experimentally relevant case 
of stretchable membranes, finding a slightly modified but 
similar form of criticality. We then confirm the existence 
of the critical scaling behavior in experiments and uti- 
lize it to obtain a reliable measurement of permeability 
coefficients. 



A crossover between two stages of vesicle swellingis 
well known in micropipette-aspiration experiments [8[. 
The stages are distinguished by a markedly different de- 
pendence of surface tension on surface strain — in the 
first stage the dependence is exponential, while in the sec- 
ond it is linear [cHlQf. The essential difference between 
this scenario and the one addressed here lies in the differ- 
ent control parameters. In the former case the swelling is 
controlled by a hydrostatic pressure difference, whereas 
in the latter it is controlled by the number of encapsu- 
lated solute molecules. As a result, the two crossovers are 
not equivalent. As will be presented below, in fact, the 
monitored swelling of our vesicles, with its two distinct 
stages, occurs within the linear regime — i.e., the ten- 
sion depends linearly on strain throughout the observed 
transition. 



We begin by describing the experimental setup in Sec. 
ILTl and the direct experimental results in Sec. IIIII The 
analysis of the experimental data requires a revised the- 
ory, which is presented in Sec. IIV1 In Sec. [V] we apply 
the theory to the experimental results to demonstrate 
data collapse onto a universal curve — i.e., the law of cor- 
responding states — and extract additional information, 
such as the permeability coefficients and the number of 
surface bending modes contributing to vesicle thermody- 
namics. Finally, we discuss the various results and their 
significance in Sec. IVII 
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II. EXPERIMENTAL SETUP 

D-(+)-glucose, D-(+)-sucrose, glycerol, urea, and 
ethylene glycol were purchased from Fluka (Buchs, 
Switzerland). Methanol and chloroform were purchased 
from Kemika (Zagreb, Croatia), l-palmitoyl-2-oleoyl-sn- 
glycero-3-phosphocholine (POPC) was purchased from 
Avanti Polar Lipids (Alabaster, USA). All the solutions 
were prepared in double-distilled sterile water. 

A suspension of POPC giant unilamellar vesicles 
(GUVs) in 0.1 or 0.2 mol/L 1:1 sucrose/glucose solu- 
tion was prepared using an electroformation method, de- 
scribed in Ref . ITT1 with some modifications [H|,[l3j]. Lipids 
were dissolved in a mixture of chloroform/methanol (2:1, 
v/v) to a concentration of 1 mg/mL. A volume of 25 fiL 
of the lipid solution was spread onto a pair of Pt elec- 
trodes and dried under reduced pressure (water aspira- 
tor; w 60 mmHg) for 2 hours. The electrodes were then 
placed into an electroformation chamber, which was filled 
with 0.1 or 0.2 mol/L sucrose. AC current (8 V, 10 Hz) 
was applied, and the voltage and frequency were reduced 
in steps to the final values of 1 V and 1 Hz 13]. Sub- 
sequently, the chamber was first drained into a beaker 
and then flushed with an equal volume of isomolar glu- 
cose solution, thus resulting in a suspension of GUVs 
containing entrapped sucrose in a 1:1 sucrose/glucose so- 
lution, which increases the contrast in a pha se contrast 
setup and facilitates vesicle manipulation |14j . This pro- 
cedure yields mostly spherical unilamellar vesicles, with 
diameters of up to 100 /jm. 

An inverted optical microscope (Nikon Diaphot 200, 
objective 20/0.40 Ph2 DL) with micro-manipulating 
equipment (Narishige MMN-l/MMO-202) and a cooled 
CCD camera (Hamamatsu ORCA-ER; C4742-95- 
12ERG), connected via an IEEE-1394 interface to a 
PC running Hamamatsu Wasabi software, was used to 
obtain phase contrast micrographs. In the streaming 
mode, the camera provides 1344 x 1024 12-bit grayscale 
images at a rate of 8.9 images/s. 

In the experiment, an individual spherical GUV is se- 
lected, fully aspirated into a glass micropipette whose 
inner diameter exceeds the vesicle's diameter, and trans- 
ferred from a solution containing solutes of very low 
membrane permeability (1:1 glucose/ sucrose) into an iso- 
osmolar solution of a more permeable solute (glycerol, 
urea, or ethylene glycol), where the content of the mi- 
cropipette is released, and the micropipette is subse- 
quently removed. Vesicle response is recorded using a 
CCD camera mounted on the microscope. The radius 
of the vesicle's cross-section was determined from the 
recorded series of images using a least-squares procedure 



III. EXPERIMENTAL RESULTS 

A transient increase of the vesicle's cross-sectional ra- 
dius, is observed immediately after its transfer. We 



attribute it to a small hypotonicity of the outer solution, 
causing the vesicle to slightly deflate and change its shape 
from a sphere into an oblate spheroid 0, [lTj • Gravity 
should lead to a small deviation of the shape from a per- 
fect oblate spheroid 16], which we neglect here. A more 
important effect of gravity is the breaking of the prob- 
lem's rotational symmetry, such that the observed lateral 
radius, R\, is always the spheroid's larger radius. On a 
longer time scale than the transient deflation, the vesicle 
inflates due to solute permeation and the accompanying 
osmosis (i.e., water influx). This causes the vesicle to be- 
come more spherical, making the observed R\ decrease. 
(See Fig.[TJ) This is the ironing stage, where the swelling 
increases the volume-to-area ratio. After a certain time, 
R\ reaches a minimum and starts increasing. (See again 
Fig. lU) This marks the crossover to the stretching stage, 
where the vesicle continues to swell essentially as an in- 
flating sphere. When the membrane reaches a critical 
strain, it ruptures, and the vesicle bursts. Subsequently, 
the membrane is resealed and another cycle of swelling 
commences. The repeated burst-swelling cycles were an- 
alyzed elsewhere [13| . In the present work we focus on 
the swelling that precedes the first burst. 

The ironing and stretching stages are clearly distin- 
guished in Fig. [1] as the decreasing and increasing parts 
of the swelling curve, respectively. To have a systematic 
definition of the time axis for all vesicles we define t = 
at the minimum of the curve. We associate the radius at 
that minimum, Rq, with the relaxed area of the vesicle, 
Aq = 4ttRq. (We return to examine this assumption in 
Sec. EE ) From each such curve we also extract the ra- 
dius at rupture, i? rup . This yields a direct measurement 
of the membrane's rupture strain, 

Crup = {Rrup/Ro} 2 — 1- (1) 

Figures HKa)-(c) show the swelling curves measured 
for POPC vesicles of broadly distributed sizes due to 
the permeation of three different solutes: urea, glycerol, 
and ethylene glycol. Each of these curves has the typi- 
cal shape shown in Fig. [TJ yet the polydispersity of the 
vesicles flattens them once they are displayed together. 
Detailed analysis of the temporal swelling curves requires 
a theory, which is presented in the next section. In Fig. 
[5] it is apparent, nonetheless, that there is no correlation 
between the transition width (sharpness of the curve's 
minimum) and the vesicle size. We return to this point 
in Sec. El 

A dataset that can be presented without further anal- 
ysis is the one for e rup , the critical surface strain at rup- 
ture, calculated from the directly measured values of Rq 
and i? rup according to Eq. (JXJ) . Since rupture is a result 
of pore formation, which is a nucleated event [III fl9| . 
the critical strain as measured in individual vesicles is 
a stochastic variable. Figure [3] shows the histogram of 
measured rupture strains. From this distribution we get 
£rup = 0.055 ± 0.02, i.e., about 5 percent, which agrees 
with known values for the maximum sustainable strain 
of lipid bilayers [|[ . 
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FIG. 1. Typical swelling curve, showing the change in the 
observed (larger) radius of an oblate spheroidal POPC vesicle 
as a function of time due to inward permeation of urea. The 
decreasing part of the curve corresponds to the ironing stage, 
and the increasing one to the stretching stage. The minimum 
of the curve is defined as (t = 0, Ri — Ro)- The present 
analysis concerns the swelling up to the point of first rupture 

(Rl — -Rrup)- 



IV. THEORY 

Models of vesicles usually focus on the statistical me- 
chanics of the membrane under the appropriate con- 
straints, the star ting point being the Helfrich Hamilto- 
nian [l|, 0, H, EH Ell]. I n the current work we choose a 
different, more general description, based on the ther- 
modynamics of small systems. This is because the phe- 
nomenon of main interest here — the transition between 
the ironing and stretching stages — arises from competi- 
tion between volume and surface effects, which is inde- 
pendent of the detailed interactions. Since the transition 
is of a mean-field type Q , such a thermodynamic descrip- 
tion should be sufficient. 

We consider a dilute solution enclosed in a membrane 
vesicle. This system is in thermal and pressure contact 
with a much larger solution having temperature T and 
pressure p . The system is at thermodynamic equilib- 
rium. Thus, its temperature is T . However, since we 
deal with strong surface effects, the system's equilibrium 
pressure may deviate from the outer pressure p . The en- 
closed solution contains Q solute molecules Q Although 
Q in the experiment changes with time, we assume, be- 
cause the solute has much lower permeability than water, 
that the system is at all times in quasi-equilibrium, given 
the instantaneous number of enclosed solute molecules, 
Q(t). Due to the fast water exchange, we treat both the 
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FIG. 2. Swelling curves for POPC vesicles under the perme- 
ation of (a) urea, (b) glycerol, and (c) ethylene glycol. Outer 
solute concentrations are 0.1 M, 0.2 M (left-pointing arrows), 
and 0.11 M (right-pointing arrow in (b)). 



In the experiment described above the vesicle actually contains 
two solutes — one that permeates the membrane and another 
that is trapped inside the vesicle. Analysis for an ideal solution, 
containing Q\ molecules of the former and Q2 molecules of the 
latter, leads to identical results once we take Q = Qi + Q2- 



volume V and surface area A of the vesicle as free, in- 
dependent thermodynamic variables. Alternatively, one 
can define a reduced volume, 

u (y,^) = 6^F 1 Le(0,i), (2) 
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FIG. 3. Histogram of measured values of rupture strains for 
POPC vesicles. The distribution has a mean of 0.055 and 
standard deviation of 0.02. 



which characterizes the deviation of the vesicle from a 
perfect spherical shape, and consider v and A as the ther- 
modynamic variables. 

We divide the Gibbs free energy of the system into 
three-dimensional (volume) and two-dimensional (sur- 
face) contributions: 



G = G 



3D 



2D- 



(3) 



Taking the enclosed solution to be ideal and dilute, we 
write its Gibbs free energy as 



G 



3D 



k B TQ 



In 



V 



1 



■PoV, 



(4) 



where v w is the molecular volume of water. 

As to the surface, we distinguish between in-plane 
(stretching) and out-of-plane (bending) degrees of free- 
dom, and assume that the vesicle is sufficiently large such 
that the two are decoupled, i.e., contribute additively to 
the Helmholtz free energy, 



F2D = F s + Fb. 



(5) 



We further assume that the membrane's in-plane free en- 
ergy is minimum for a certain relaxed area, Aq, and ex- 
pand it to leading order about that minimum, 



A,) 2 , 



(6) 



where K is the membrane's stretching modulus. In the 
bending free energy we include only the contribution from 
undulation entropy, which, assuming small fluctuations 
about a spherical shape, is given by [7| 



F h = -^Nk B Tki(l-v), 



(7) 



where N is the total number of independent bending 
modes contributing to the membrane's thermodynamics. 



Equation (J7J, which describes the suppression of bend- 
ing fluctuations as the vesicle approaches the maximum 
volume-to-area ratio (v — > 1), is the only statistical- 
mechanical input to the theory. It is analogous to an 
equipartition principle; here each bending mode con- 
tributes an equal amount of (k-BT/2)\ ln(l — v)\ to the 
free energy [7|. Thus, this expression holds quite gener- 
ally for nearly spherical vesicles, regardless of their de- 
tailed properties. Equation ([7J also provides the crucial 
coupling between G3D and F2D , since v depends on both 

V and A according to Eq. ([2]). 

The number of modes affecting the thermodynamics, 
N, which should be proportional to the number of lipid 
molecules in the membrane (i.e., to the membrane's re- 
laxed area Aq), is a priori unknown. Further discussion of 
this number, its dependence on membrane's parameters, 
and its link to other theoretical treatments of membrane 
fluctuations, is deferred to Sec. IVD Here we prefer to 
leave it unspecified and extract it from the experiment. 
The number of modes is used to define an intensive area, 
or effective patch size, a = A/N, having the relaxed value 
a Q = A /N. 

One could add, of course, various other contributions 
to the free energy of the vesicle, such as the membrane's 
bending energy, surface-surface and surface-solute in- 
teractions (e.g., electrostatic ones), gravitational energy, 
etc.. All these, however, do not significantly change when 
the vesicle becomes spherical and begins to stretch. In 
other words, in the limit (to be studied below) where 
(1 — v) becomes singular, these contributions are non- 
singular, and, hence, cannot affect the critical behavior 
in the transition region. 

Equations (J2])-© define the Gibbs free energy of our 
finite-size system as a function of (T,p Q ,Q, N) and also 

V and A. The equilibrium free energy as a function of 
(T,p , Q, N) alone is obtained by minimizing G with re- 
spect to the volume and area of the vesicle. Note that by 
using this procedure we circumvent altogether the sub- 
tle and controversial issue of membrane tension. (See, 
e.g., Refs. I20U22I and references therein.) The Laplace 
tension and surface pressure of the membrane emerge 
naturally from the minimization with respect to V and 
A, respectively. That is also the reason why we have for- 
mulated the Helmholtz (rather than Gibbs) free energy 
of the surface, to avoid specifying any thermodynamic 
variables pertaining to the surface other than T, A, and 
N. Another relatd subtlety is that K = A(d 2 F 2D /dA 2 ) 
at fixed v (rather than fixed V), otherwise there would 
be another contribution to the stretching modulus from 
changes in the volume-to-area ratio. 

The minimization leads to the following two equations: 
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(8) 

(9) 



from which one can calculate v and a. In Eqs. © and 
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© we have defined 



PqVq 
k B T' 



A 



3/2 



along with two small parameters, 
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(10) 



(11) 



the first related to the vesicle's finite size, and the other 
to the membrane's finite stretchability. 

Before proceeding, it is instructive to examine Eqs. 
© and © in two limits. The first is the ordinary 
thermodynamic limit of infinite system size (N — > oo, 
<5jv — > 0) while keeping the stretching modulus finite. In 
this case Eqs. © and © become degenerate and yield 
the expected equilibration of inner and outer pressures, 
k B TQ/V — p a . The second limit is that of an unstretch- 
able membrane (K — > oo, 6k — > 0). In this case Eq. © 
imposes a — ao, and Eq. then yields the behavior that 
was studied in Ref. 0, with criticality of v as N — » oo. 
Outside these two limits we get from these equations the 
relation 



•K 



4a(a — ap) 



3a§ 



(12) 



which demonstrates the actual interplay between the 
volume-to-area ratio and surface strain. We see that the 
vesicle can never attain a perfect spherical shape (v = 1), 
because this would require unphysical strain (a — > oo). 
The deviation of the area from its relaxed value is in- 
versely proportional to the deviation of the shape from 
a perfect sphere. This relation becomes more sensitive 
the smaller the value of 6 k — i.e., the vesicle should be 
strongly swollen, 1 — v ~ Sk "C 1, to get an appreciable 
strain (a — asoV&o- 

To present the detailed behavior of the transition, we 
define a control parameter q, proportional to the number 
of enclosed solute molecules, and an order parameter M, 
measuring the deviation from a spherical shape, 



q=Q/Q c -l, M=l 



(13) 



Solving Eqs. 
we get 



and (JSJ) in the vicinity of the transition, 



M{q) = AM (g/A) 

m{x) = (yr 



A = [26 N 
x) /2. 



(9/2)6 



K 



nl/2 



(14) 



Thus, as q increases from negative to positive values, M 
crosses over from appreciable values, M ~ \q\, to very 
small ones, M ~ A/(4g). The crossover occurs over a 
range of q defined by A. In the limit A — > the point 
(q = 0, M = 0) becomes a singular corner where dM/dq 
undergoes a discontinuous jump. The finite size of the 
vesicle and the stretchability of the membrane both con- 
tribute to the broadening of the transition. If we take 



6k = 0, Eq. (TT4"|) properly coincides with the critical 
behavior in the case of an unstretchable membrane 
Which of the two broadening effects dominates in prac- 
tice is determined by the ratio 



6n 
Sk 



K/R 



Po 



where the right-hand side reflects the competition be- 
tween surface (numerator) and volume (denominator) ef- 
fects. 

Substituting Eq. ([14]) in Eq. ([12]), we obtain for the 
surface strain, e = (a — a )/a Q = (3/4)<5^M _1 . Yet, 
since the scaling function satisfies [M (a;)] -1 = 4M(— x), 
we can rewrite this result as 

- wmw^A Am - q/A) - (15) 

Thus, up to a prefactor, the strain's behavior is a mirror 
image of that of M. As q increases from negative to pos- 
itive values, e changes from very small values, ~ A/|q|, 
to appreciable ones, ~ q. Similarly, combining Eqs. ([8]), 
([T2"]). and ([Tf]). we get for the Laplace tension, 



7(9) 



R 
~2 



k B TQ 
V 



Po 



3k B T 
4a n 



[M{q)]' x = Ke(q), 



i.e., the surface tension follows the same behavior as the 
strain, becoming appreciable only for sufficiently large, 
positive q. Equation (fT6)) implies that within the tran- 
sition region the tension is in the so-called linear regime 
[H Q , where it is proportional to the surface strain. Sub- 
stituting the rupture strain found above, e rup ~ 0.05, and 
K ~ 240 mN/m [23|, we obtain a rupture tension of or- 
der 10 mN/m, which is in line with previously reported 
values 

Equation (fl"4]) [equivalently, (fT5]l or (IT6l) ] describes 
a law of corresponding states for osmotically swollen, 
nearly spherical vesicles. It predicts that, upon rescal- 
ing, the swelling (i.e., volume-to-area ratio, strain, or 
tension) of such vesicles, regardless of their size, compo- 
sition, and the nature or concentration of solute, could 
be collapsed onto a single universal curve defined by the 
function M(x). 



V. DATA ANALYSIS 

Let us first check which of the broadening factors — 
finite size or stretchability — is dominant in the exper- 
iment. The stretching modulus of a POPC bilayer is 
K ~ 240 mN/m . Our vesicles have R - 20-50 /xm, 
and the outer concentration is c ~ 0.1-0.2 M. We get 
6 N /6 K - K/(R k B Tc o ) - 0.01-0.1. Hence, the transi- 
tion width is governed by the finite stretching modulus 
and is expected, therefore, to be independent of vesicle 
size, 



A 



-8k 



1/2 



= ( 9k B T \ 
\2KaJ 



1/2 



(17) 
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To further check the theoretical predictions of the pre- 
ceding section, we need to relate them to the experimen- 
tal observables demonstrated in Fig.[T] — i.e., the tempo- 
ral change in the principal radius of an oblate spheroidal 
vesicle, Ri(t). 

We begin with the time axis. The membrane perme- 
ability coefficients of the examined solutes are of order 
P ~ 1CT 2 fim/s H 

EH HII' the vesicles are of radii 
Ro ~ 20-50 fim, and the swelling process lasts about 
t ~ 10-100 s. We have (P/R )t ~ 0.001-0.05, implying 
that the concentration of permeating solute inside the 
vesicle never exceeds a few percent of its outer concen- 
tration. Hence, throughout the observed swelling we may 
assume dQ/dt = PAqc q , leading to 



q = (3P/R )t + const. 



(18) 



Thus, the time axis, up to a linear transformation that 
depends on vesicle size, is equivalent to our control pa- 
rameter. 

Next, we should adapt the theoretical scaling relations, 
Eqs. (TH|) and (fTo)) . so that they could be applied to the 
measured swelling curves, Ri(t). This calculation is pre- 
sented in the Appendix. Given the experimentally mea- 
sured dependence, [Ri(t)/R ] 3 — 1 = g C x P (t), we find that 
the following scale and shift transformation: 

/(*) = (15A/8)- 1 / 2 {g exp [x - (15^ 3 /A)A-^ 3 } - xA 
+3(15 1/3 /4)A 2 / 3 } , (19) 

where x = q/A and q = (3P/Ro)t, should collapse the 
data onto the universal function 



fix) = Wl + x^ -x) 1 ' 2 . 



(20) 



The rescaling scheme defined in Eq. (TT9")) , which has 
been dictated by the nature of the experiment, is not as 
elegant as the theoretical law of corresponding states, Eq. 
([T4|). Nevertheless, the two are equivalent and similarly 
straightforward. We have applied the scheme to the data 
of Fig. [5] while using the permeability P (which merely 
scales the horizontal axis) and the transition width A as 
two fitting parameters. The results are presented in Fig. 
IH showing successful data collapse, within the transition 
region and above it, onto the predicted universal curve, 
Eq. (HI]). 

In Fig. [5] we present the fitted values for the POPC- 
membrane permeability coefficients of the three exam- 
ined solutes. The values for the permeability of urea 
[Fig. IH a )] are narrowly distributed and concentration- 
independent, yielding P = 0.013 ± 0.001 ^m/s. This 
agrees well with the result of 0.014 fim/s, obtained for 
0.1 /im-radius DOPC vesicles by Paula et al. 5] using 
dynamic light scattering, despite the different phospho- 
lipid and the large difference in the sizes of the stud- 
ied systems. For glycerol we get values that signifi- 
cantly increase with glycerol concentration [Fig. [5jb)]: 
P = 0.0053, 0.0074, and 0.019 ± 0.006 /Ltm/s for c Q = 0.1, 



(a) 

Q 
o 


1 

\ 


1 

urea 


6 








~ X - \v\\ 




2 



- ^ 




1,1,1,1, 



-40 



-20 



20 



40 



60 



(b) 



X 



glycerol 




-10 



10 20 
X 



j i z 

30 40 



(C) 



X 



ethylene glycol 




-20 



FIG. 4. The experimental data of Fig. [2] replotted using the 
same colors after being transformed according to Eq. 
The dashed gray line shows the theoretical universal curve, 
f(x) = (VT+^-x) 1/2 . 



0.11, and 0.2 M, respectively. The last value agrees with 
the value of 0.021 ± 0.008 (im/s, extracted for the same 
vesicles, at c = 0.2 M glycerol, from a different analysis 
of the stretch-burst cycles that follow vesicle rupture [l3| . 
A permeability of 0.027 /im/s was reported by Paula et 
al. [5J and by Dordas and Brown [24| for 0.1 /zm-radius 
egg-PC vesicles at c Q = 0.4 M glycerol, which is not in- 
consistent with our results, especially given the concen- 
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tration dependence. For ethylene glycol we obtain more 
broadly distributed values that depend on concentration 
[Fig. BUM]: P = 0.046 ±0.006 and 0.085 ±0.01 ^m/s for 
c = 0.1 and 0.2 M, respectively. The only previously 
measured permeability for ethylene glycol that we are 
aware of is 0.88 /im/s, obtained by Orbach and Finkel- 
stein [25j using conductivity measurements for a flat egg- 
PC membrane. This value is larger than ours by an or- 
der of magnitude. (We note, though, the very different 
conditions under which the two experiments were con- 
ducted.) All permeability values reported here, includ- 
ing the concentration-dependent ones for glycerol and 
ethylene glycol, are consistent with those found from the 
alternative analysis of stretch-burst cycles (unpublished 
results) . 
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FIG. 5. Permeability coefficients as a function of outer solute 
concentration for urea (a) and all three solutes (b). The lines 
are guides to the eye. 



tic variable. Of the three parameters appearing in Eq. 
(TT7)l — T, K, and an — only the effective patch size an 
may be responsible for such stochasticity. In Fig. |6fb) 
we show the distribution of patch sizes, as arising from 
Eq. (JT7J) using the measured A and K = 240 mN/m. The 
mean patch size is an ~ 0.01 /^m 2 , which is much smaller 
than the vesicle area, but much larger than the molecu- 
lar size. Such a patch contains about 10 4 -10 5 lipids. In 
the inset of Fig. Mb) we recast the same data in terms of 
the number of effective bending modes, N = An/an. The 
number of modes is broadly distributed in the range 10 5 - 
10 7 , making its increasing trend with vesicle area hardly 
discernible. 
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FIG. 6. (a) Transition width vs. vesicle radius, exhibiting 
no correlation between the two variables, (b) Histogram of 
effective patch size values, obtained from the measured A of 
panel (a) using Eq. (|17[) and K = 240 mN/m. Inset shows the 
resulting number of effective bending modes vs. vesicle area. 



We now turn to the second fitting parameter — the 
transition width, A. It is found to lie in the range 10 -3 - 
10~ 2 and be uncorrelated with the vesicle size (Fig. [5]), as 
predicted above. The small values of A establish the self- 
consistency of our analysis for the crossover between the 
two swelling stages as a slightly broadened phase transi- 
tion, and thus explain the success of the resulting data 
collapse. The rather broad distribution of these values 
imply that the transition width is affected by a stochas- 



VI. DISCUSSION 



The experimental results presented above confirm the 
critical nature of the osmotic swelling of vesicles, as was 
first suggested in Ref. 0- It is important to note that 
a similar critical scaling is absent in the case of hydro- 
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static swelling, as in vesicle aspiration experiments [8j. 
To adapt our formulation to such a case we should merely 
replace the solute contribution to the Gibbs free energy 
(the first term in Eq. Q), with (—PmV), where p- ln is the 
inner hydrostatic pressure. Repeating the minimization 
procedure presented in Sec. IIV1 we find that Eq. (fT2|) 
remains unchanged, while the swelling order parameter 
becomes 

M = l-V=- / 7^ -r. (21) 

2(j3 in / Po -l) v ' 

Thus, our formalism applies in this case only for suffi- 
ciently large inner pressures, Pi n /p > I+Sn/2 (to ensure 
a positive v). Similar to q = Q/Q c — 1, we may define a 
swelling control parameter, q' — p ln /p —l 7 yet its value is 
restricted to q' > Sn/2. Consequently, Eq. (j2Tj) yields a 
moderate decrease of M from appreciable to small values, 
M ~ 1/q' — i.e., a gradual approach to a spherical shape. 
Another way to demonstrate the non-criticality of the 
hydrostatic-pressure case is to recognize that, once both 
p and p m are specified, the surface tension of a nearly 
spherical vesicle is simply given by Laplace's law rather 
than the more elaborate critical expressions of Eqs. (|14D - 
(|16[) . This underlines the special thermodynamic behav- 
ior of osmotically stressed vesicles, as was indicated be- 
fore in a more general statistical-mechanical context (2(| . 
The key difference between osmotically swollen and di- 
rectly pressurized vesicles lies in the fact that the pres- 
sure difference of the former depends on the volume for 
a given number of encapsulated solute molecules [cf. Eq. 
(UJ)], whereas the pressure difference of the latter is a 
given constraint. This makes the statistical ensembles of 
(T,p ,Q,N) and (T,p ,p; n , N) not equivalent [26j . 

We have shown how the data collapse for different vesi- 
cle sizes can be utilized to obtain a sensitive measurement 
of the membrane permeability coefficients of various so- 
lutes. This ability is particularly well demonstrated in 
the case of urea (Fig. Ufa)). Overall, the permeability 
values that we have obtained are in accord with previ- 
ously published ones. (See the detailed account in the 
preceding section.) In addition, we have found an in- 
creasing concentration dependence of the permeability 
for glycerol and ethylene glycol. While we are not aware 
of an earlier report of such an effect, it is in line with 
strong evidence for the affinity of polyols to lipid head- 
groups [13, and may account for the large range of 
values reported in the literature for the permeability coef- 
ficients of these solutes. The affinity should make the so- 
lute concentration adjacent to the membrane larger than 
in the bulk solution, leading to an apparent permeabil- 
ity coefficient which is larger than the actual one. More 
extensive measurements of the apparent permeability as 
a function of concentration might actually allow an ex- 
traction of the solute-membrane affinity parameters. At 
sufficiently high concentrations the polyols may also dis- 
rupt the membrane structure and thus affect the actual 
permeability p9| . 

The transition width, and the corresponding effective 



membrane patch oo, have been found to be broadly dis- 
tributed (Fig. [5]). The origin of this stochasticity is not 
yet clear. It may be that ao is sensitive to the pres- 
ence of membrane impurities. Nevertheless, the analysis 
of transition widths has enabled us to obtain information 
regarding the number of independent bending modes, N, 
which contribute to the thermodynamics of the vesicle in 
the transition region. To our knowledge this variable has 
not been experimentally accessed before. For the vesi- 
cles used here, having Rq ~ 20-50 /zm, we have found 
numbers N ~ 10 5 -10 7 , which are large but much smaller 
than the number of lipid molecules in the vesicle (about 
10 10 ). 

To account for these intermediate values of N we 
should link our thermodynamic formulation to the more 
detailed statistical mechanics of fluctuating vesicles. The 
definition of N in Sec. IIVI has been through the sur- 
face entropy of a nearly spherical vesicle, Eq. ([7]), where 
each of the N modes contributes an equal amount of 
(kBT/2)\ ln(l — v)\ to the free energy Examining 
more closely the partition function of a fluctuating vesi- 
cle within Helfrich's model [l(|, we find that these modes, 
whose contribution is singular in (1 — v), belong to the 
small-wavenumber (tension-dominated) portion of the 
fluctuation spectrum. Their number is 

N ~ jRl/n, (22) 

where k is the membrane's bending rigidity. Equations 
(fT5)l . (1P51) . and (|2"2"j) then allow a theoretical estimate of 
N. At the transition we have 7 ~ (2/3)K AM(-l) (cf. 
Fig. IU, which for K ~ 240 mN/m 23] and A ~ 0.004 
(Fig. Hl[a)) yields 7 ~ 0.8 mN/m. This value of about 1 
mN/m is in line with that obtained for nearly spherical 
vesicles in aspiration experiments @. Substituting it, 
together with k ~ 10fc B T [H and R ~ 20-50 fim, in 
Eq. , we get N ~ 10 7 , which is in the range of values 
that we have found experimentally. 

Equation (122)) implies, in fact, that N is not constant 
but increases with the degree of swelling (i.e., with 7). 
Hence, sufficiently far before the transition 7 will be too 
small, and the theory should fail. (Far beyond the tran- 
sition, on the other hand, the behavior is dominated by 
surface strain and unaffected by these entropy consid- 
erations.) Combining Eqs. (fTS|) and (|22l) . we find that 
the theory should break down for M > 3fceT 1 /(167rK) ~ 
0.006. With A ~ 0.004 this occurs for x < -1.3, which is 
consistent with the region where the data collapse begins 
to fail (see Fig. gj. 

Thus, the assumption of a swelling- independent N is 
the main limitation of the model and should be improved 
in future studies. Other assumptions are less significant. 
For example, the relaxed area of the vesicle, Aq, is actu- 
ally slightly smaller than the assumed value of 47ri?j^ in , 
where i? m i n is the minimum of the experimental R±(t) 
curve. The calculation given in the Appendix implies 
that Ro/Rmin ~ 1 - (15 1 / 3 / 4 ) A2/3 - For OU1 ' broadest 
transitions this deviation amounts to less than 3%. As 
shown in Ref . 0, accounting for additional factors such as 
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a non-ideal solute does not affect the universal behavior 
in the vicinity of the transition. 
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APPENDIX 

In this appendix we adapt the theoretical scaling rela- 
tion obtained in Sec. HVI to the observables in the experi- 
mental swelling curves. The result is the scaling transfor- 
mation that should be applied to the experimental data 
to achieve data collapse. 

Consider an oblate spheroid of principal radii R\ and 
i? 2 < Ri, and small eccentricity e = [1 - (i? 2 / Ri) 2 ] 1 ^ 2 <C 
1. The volume of the spheroid is V = {A-k/Z)R\ R 2 , 
and its area is A = it{2R 2 + (i2§/e)ln[(l + e)/(l - e)]}. 
From this we get v ~ 1 — e 4 /15 + 0(e 6 ), leading to 
e 2 ~ (15A/) 1 / 2 + O(M), and R 2 /R x = (1 - e 2 ) 1 / 2 ~ 
l-(15M) x /2/2 + 0(M). Now we have, on the one hand, 
V/V = (A/A ) 3 / 2 v = (1 + e) 3/2 v, and, on the other 
hand, V/Vq = {Ri / R ) 3 {R 2 / Ri) ■ Equating these two 
expressions while using the result above for R 2 /Ri, we 
obtain 

{Ri/Rof - 1 ~ (15M) 1/2 /2 + 3e/2. (Al) 

The correction to this expression is 0(M). In the tran- 
sition region this correction, as well as e, are negligible 



(~ A) compared to the M 1 / 2 term (~ A 1 / 2 ). Ab ove the 
transition, at the point where e becomes comparable to 
M 1 / 2 , we find e - M 1 ' 2 ~ A 2 / 3 , whereas M ~ A 4 / 3 . 
Thus, Eq. (|Aip contains the leading terms within the 
transition region and above it. Using Eqs. (fT4")) , (|15p . 
and the property M(— x) = M(x) + x, we rewrite Eq. 
dAl~|) as 

(R 1 /R ) 3 -l~(15A/8) 1 / 2 f(q/A) + q 

f(x) = (^T+x^-x)^ 2 , (A2) 

where we have omitted another term of order M and also 
the ratio 8n/v~k- 

Equation (|A2I) . as required, contains the experimen- 
tal observable R\/Rq rather than our order parameter 
M. Yet, we need to take one last step before we can ap- 
ply it to the experimental data. The expression in Eq. 
([MJ) has a minimum at [q ~ (15 1 / 3 /4)A 2 / 3 , (i?i/i? ) 3 - 
1 ~ 3(15 1/3 /4)A 2 / 3 ], which is shifted from the one 
that we have defined for the experimental curve, [q = 
0, (R1/R0) 3 -1 = 0] (cf. Fig. EE). In other words, the 
actual relaxed radius is slightly shifted from the mini- 
mum of the measured swelling curve. To correct for this 
~ A 2 / 3 error, we ought to shift the experimental curves 
accordingly. If [Ri(t) j Ro} 3 — 1 = g C x P (t) is the measured 
dependence, as a function of q, then the following scaling 
and shifting transformation, 

f(x) = (15A/8)- 1 / 2 { 3cxp [x - (15 1 / 3 /4)A- 1 / 3 ] - xA 

+ 3(15 1 / 3 /4)A 2 / 3 } 

x = q/A, q = {3P/Ro)t, (A3) 

should collapse the data onto the universal function f(x) 
defined in Eq. (IA2[) . This is the transformation used in 
Sec. El 
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